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The quasiparticle states around a nonmagnetic impurity in electron-doped iron-based superconductors with 
spin-density-wave (SDW) order are investigated as a function of doping and impurity scattering strength. In the 
undoped sample, where a pure SDW state exists, two impurity-induced resonance peaks are observed around the 
impurity site and they are shifted to higher (lower) energies as the strength of the positive (negative) scattering 
potential (SP) is increased. For the doped samples where the SDW order and the superconducting order coexist, 
the main feature is the existence of sharp in-gap resonance peaks whose positions and intensity depend on 
the strength of the SP and the doping concentration. In all cases, the local density of states exhibits clear C2 
symmetry. We also note that in the doped cases, the impurity will divide the system into two sublattices with 
distinct values of magnetic order. Here we use the band structure of a two-orbital model, which considers the 
asymmetry of the As atoms above and below the Fe-Fe plane. This model is suitable to study the properties 
of the surface layers in the iron-pnictides and should be more appropriate to describe the scanning tunneling 
microscopy experiments. 

PACS numbers: 74.70.Xa, 74.55.+V 



I. INTRODUCTION 

The new family of iron-based superconducting (SC) mate- 
rials has attracted much attention since their discovery.[3 The 
parent compounds exhibit a spin-density-wave (SDW) order 
at low temperatures. Upon doping either electrons or holes 
into the system, the SDW order is suppressed and supercon- 
ductivity emerges, suggesting the interplay and competition 
between these two states. 

The impurity effect is an important property in the studies 
of superconductivity. One prominent feature of li-wave par- 
ing symmetry in cuprates is the existence of bound states near 
the impurity, which is revealed by both experiments and the- 
oretical calculations.- For iron-based materials, the impurity 
effect in the SC state has also been theoretically studied inten- 
sively.^' '° It was proposed that a single nonmagnetic impurity 
could be used to distinguish the pairing symmetry and the in- 
gapbound states could exist for the typical s±- pairing symme- 
tryPuni'nie in-gap bound states for Sj- pairing should be differ- 
ent from those in cuprates due to the absence of quasiparticle 
excitations at low energy. On the other hand, in the SDW state, 
it was proposed experimentally that the Fermi surface (FS) is 
only partially gapped and small ungapped Fermi pockets exist 
at low temperature.'^^EIIxhis feature was recently reproduced 
based on a two-orbital model together with a mean-field ap- 
proach.'^ The gap-like feature and the existence of tiny un- 
gapped regions along the diagonal direction of the Brillouin 
zone (BZ) are quite similar to the case of the li-wave SC gap 
in cuprates. Therefore, one would expect that the low-energy 
bound states should also exist in the SDW state for the iron- 
based materials. So far, the impurity effect in the SDW state 
remains less explored and a systematic study for this topic is 
still lacking. We believe it is timely and quite interesting to 
address this issue theoretically and verify the aforementioned 



expectation numerically. Furthermore, in some of the iron- 
based materials the SC and SDW orders are proposed to coex- 
ist in the underdoped regime.'^^'^ The issue is still a subject 
of discussion and we anticipate that the impurity effect could 
provide additional signatures for the coexistence of these two 
orders. 

In this paper, we study theoretically the impurity effect on a 
two-dimensional square lattice based on a two-orbital model 
and the Bogoliubov-de Gennes (BdG) equations. By introduc- 
ing a single impurity into the system, the local density of states 
(LDOS) is calculated and our results show that: (i) In the 
undoped sample, there are two impurity-induced resonance 
peaks at and near the impurity site and the LDOS spectra ex- 
hibit C2 symmetry, with one-dimensional modulation, (ii) The 
impurity effect in various doped cases are also studied. Its ef- 
fect on the LDOS is remarkable only when the strength of the 
scattering potential (SP) is larger than a certain value. For 
weak and moderate SPs, a distinct bound state exists explic- 
itly at the next-nearest-neighbor (nnn) sites of the impurity, 
the energy of which depends on the strength and sign of the 
SP, as well as on the doping concentration. For the unitary 
impurity, there is a sharp in-gap peak at low doping; while at 
high doping, the impurity induced bound state is close to the 
SC coherence peaks. On the other hand, in a small range of 
moderate doping there are two in-gap peaks only for positive 
SP. All the above features could be used to detect the presence 
of the SDW order as well as the coexistence of the SDW and 
SC orders. 

The paper is organized as follows. In Secjllj we introduce 
the model and work out the formalism. In Sec. Ill we study 



the FS. In Sec. IV the impurity effect in the parent compound 



is investigated. In Sees. [V] and VI we study the impurity 
effect in doped regime for positive and negative SPs, respec- 
tively. Finally, we give a summary in Sec. VII 
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II. MODEL AND FORMALISM 



The iron-based superconducting materials have a layered 
structure with the FeAs layers being the superconducting 
planes. In the present work, following Refs. [9| and |17| 
and taking into account the asymmetry of the As ions located 
above and below the Fe-Fe plane, we start from a two-orbital 
model with on-site interaction. This model is able to qualita- 
tively explain both the ARPES^^ and STM"** experiments. For 
example, based on this model, the obtained phase diagram,'^ 
spin susceptib ility ,E2l as well as the Andreev bound state inside 
the vortex core^are all consistent with the experiments. Thus 
we also use it to investigate the impurity effect. The Hamilto- 
nian is written as. 



(1) 



Here Hgcs is the BCS-like Hamiltonian, which includes the 
hopping term and the pairing term, expressed by, 

HbCS - - Tjitiiva-itifiivcl^ijCjvcr + h.C.) - t{) ^i/jo" '^i^tr'^i/'iT 
+ 2iMV^r(Ai;,jvc]^^c]"v<r + ' 

where i - {ixJy), j - ijx,jy) are the site indices, yu, v - 1,2 
are the orbital indices, and to is the chemical potential. //,„, is 
the on-site interaction term. At the mean-field level, it can be 
written asPEUll 



Hint = ^ X ^"i/^*^"'/"^ + ^' ^'^vs-)nivo- 

(3) 



where nxf^a- is the density operator at site i and orbital ^, with 
spin cr. The quantity U' is taken to be U-2Jh-^^ The impurity 
part of the Hamiltonian, //,„,p, is given by: 



(4) 



Following Refs. II9I10I17I . we here consider only the intra- 
orbital scattering by a nonmagnetic impurity. 

The mean-field Hamiltonian (1) can be diagonalized by 
solving the BdG equations self-consistently, 



-H* . . 



u. 



(5) 



with 



Hi^iiya = -fi/ijv + Wini^a-) + (U - 2JH){nifig-} 



(6) 



and 



Here V\^jv is the pairing strength and f{x) is the Fermi-Dirac 
distribution function. The SC order parameter at site i is de- 
fined as 



A/ 



in accordance with the s± pairing symmetry. 
The LDOS is calculated according to 



x+y 



(9) 



Piicj) = Y^\ulJ^diE„ -cj) + |vV |2<J(£„ + oj)] , (10) 



where the delta function 6(x) is taken as r/7r(x^ + r^), with the 
quasiparticle damping F - 0.01. 

Following Ref. Q, we use the hopping constants. 



ti,i,i±af, = fi (a^ x,y) , 

_ 1 +(-!)' , 1 -(-!)' 

ri/i,i±(i+y)/i - 2 2 "I ^ 13 

_ 1 +(-!)' i-(-iy 

'ifi,i±{x-y)fi - 2 3 "I 2 2 

tifi,i±x±yv - h (/J. v) . 



(11) 

(12) 

(13) 
(14) 



2KbT' 



In the present work, we use fi_4 = 1,0.4,-2,0.04.121 to is 
determined by the electron filling per site n (n = 2 + x). The 
on-site Coulombic interaction U and Hund's coupling Jh are 
taken as 3.4 and 1 .3, respectively. The pairing is chosen as nnn 
intra-orbital pairing with the pairing strength V - 1.2. This 
kind of pairing is consistent with the ii-pairing^'' and has been 
widely used in previous theoretical studies based on the BdG 
technique.*^ The numerical calculation is performed on a 
32x32 square lattice with the periodic boundary conditions. A 
30 X 30 supercell is taken to calculate the LDOS. Throughout 
the paper, the energy and length are measured in units of fi 
and the Fe-Fe distance a, respectively. The temperature is set 
to be r = 0. In the following, all the results we presented 
have been checked by using different initial values and they 
remain qualitatively the same, suggesting the reliability of our 
calculation. 



III. FS IN THE SDW STATE 

In order to investigate the impurity effect in the iron-based 
superconductors with SDW order, we first study the FS in the 
SDW state. 

At zero doping, below the SDW temperature, it was pro- 
posed experimentally that there exist small FSs along the 
T-M line of the BZ and Dirac cones in the electronic structure 
form inside these FSs, with their apices being located close 
to the Fermi energy. However, whether these FSs and Dirac 
cones are electron- or hole-like is within uncertainties of the 
experiment." On the other hand, theoretically it was shown 
that in both a two-band model and a five-band model, nodes 
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FIG. 1: (Color online) (a) The zero temperature SDW FS. (b) 
Two Dirac cones at (kj,,ky) = (0.286;r, 0.286;r) and {k-,,ky) = 
(0.308;r, O.SOStt), respectively, (c) The spectral function A(k,a)) 
integrated from a> = -0.1 to toi = 0.1. (e) and (f) are the 
band structures near the Fermi energy along the blue [goes through 
(fej, ky) = (0.286;r, 0.286;r)] and orange [goes through (k^, ky) = 
(0.308;r, 0.308;r)] lines in (d), respectively. The BZ is defined in the 
2Fe/cell representation and the green dashed line in (a), (c) and (d) 
represents the MBZ. 



in the SDW gap function must exist due to the symmetry- 
enforced degeneracy at the F and M high-symmetry points, 
even in the presence of perfect nesting, but the number and 
locations of these nodes are model dependant.l^ Therefore, 
whether they correspond to the experimentally observed Dirac 
cones is still unclear. In Fig. [T] we plot the zero temper- 
ature SDW FS and the corresponding band structure near 
the Fermi energy obtained by our self-consistent calculation. 
As we can see from Fig. [T] (a), in the SDW state, there re- 
main four small FS pockets in the magnetic Brillouin zone 
(MBZ), two of which are electron-like (red) located around 
(k^, ky) - +(0.2867r, 0.2867r), while the other two are hole-like 
(blue) located around (k„ky) = +(0.3087r,0.3087r). The pock- 
ets outside the MBZ are just replica of those inside it due to 
band-folding in the SDW state and they can be connected by 
the SDW wave vector Q - (n, n). The areas enclosed by these 
pockets are equal, thus keeping the doping level at x = 0. In- 



FIG. 2: (Color online) The zero temperature SDW FS at various dop- 
ing levels. The SC order A is artificially set to zero in order to illus- 
trate the effect of SDW on the evolution of the FS. The blue and red 
pockets in the x = 0.04, 0.08 and 0. 1 cases are both electron pockets. 
The green dashed line is the same as that in Fig.[T| 



side these four pockets, there are four Dirac cones. As shown 
in Fig. [T] (b), the apex of the Dirac cone is 0.026 below the 
Fermi energy at {k^, ky) = +(0.286;r, 0.286?:) while it is 0.046 
above the Fermi energy at (A;^, ky) - +(0.3087r, 0.3087r), sug- 
gesting that they are electron- and hole-like Dirac cones, re- 
spectively. The spectral function A(k, u), which is propor- 
tional to the photoemission intensity measured in ARPES ex- 
periments, is integrated from co - -0.1 to <y = 0.1 and shown 
in Fig.[T](c). As we can see, the locations of the bright spots 
are around {k^, ky) - +(0.3?:, 0.3?:) and the equivalent symme- 
try points outside the MBZ, on the F - M line, in qualitative 
agreement with experiment. ' ' In addition, although most parts 
of the original FSs around F are gapped by the SDW order, 
the gap value is extremely small on these FSs. Thus, around 
F, the low-energy spectral function has moderate intensity and 
this can be seen from the ring structure around F with lower 
intensity, as compared to those bright spots. We also notice 
that the system has only two-fold symmetry when entering 
the SDW state while the experimentally observed four-fold 
symmetry is due to the superposition of twin domains or do- 
main averaging, as suggested in Refs. J_l,and JJ, respectively. 
The band structures near the Fermi energy scanned along the 
blue and orange cuts in Fig. [T| (d) are plotted in Figs. [T] (e) 
and [T] (f), respectively. It clearly shows the X-like struc- 
ture of Dirac cones and again suggests that the Dirac cone 
is electron-like at (A;^, ky) - {Q.2%6n, 0.286?:) and hole-like at 
{k„ky) = (0.3087r, 0.308;r). The locations of the FS pock- 
ets and the bright spots in the spectral function are consistent 
with the experimental observation, but in our calculation, the 
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electron- and hole-like Dirac cones appear in-pairs and are lo- 
cated very close to each other along the F - M line of the 
BZ, the apices of which are both in the vicinity of the Fermi 
energy, thus we propose this to be directly verified by future 
ARPES experiments with higher resolution. In addition, the 
existence of electron and hole Dirac cone states in-pairs has 
already been confirmed indirectly by measuring the magne- 
toresistance.-'^ 

Fig.[2]shows the evolution of the FS with doping. Here, we 
set the SC order A to zero to illustrate the effect of SDW on 
the evolution of the FS. In the MBZ, as doping increases, the 
size of the electron pockets [the red pockets shown in Fig. [T] 
(a)] is enlarged while that of the hole pockets [the blue pock- 
ets shown in Fig.[T|(a)] is reduced. When doping increases to 
about X = 0.02, the hole pockets vanish completely. By fur- 
ther increasing doping, another two electron pockets appear in 
the MBZ [the blue pockets in the x - 0.04, 0.08 and 0.1 cases 
shown in Fig. [2], exactly at the same locations where the hole 
pockets vanish and overlap with the original electron pockets. 
The size of all these electron pockets is enlarged with doping. 
If we define the areas enclosed by the inner and outer red lines 
to be 5 1 and those enclosed between the inner and outer blue 
lines to be 5 2, then we have x = 2N^Ny(S i + 5 2), with A^^, 
A^y being the linear dimensions of the square lattice. Finally, 
when X = 0.15, the SDW order disappears and there is no 
more band-folding due to it, in this case, there are two elec- 
tron pockets and two hole pockets around the M and F points 
of the BZ, respectively. 



IV. IMPURITY SCATTERING EFFECT IN UNDOPED 
SAMPLE 

Based on a toy model and phenomenological calculation,'^^ 
it was proposed that the impurity induced bound state should 
appear near the impurity site for the undoped sample. While 
actually the band structure and FS should be important issues 
for the features of the order parameters and LDOS. Thus we 
will reexamine this issue based on the two-orbital model and 
present a detailed investigation of the nonmagnetic impurity 
effect in the parent compound. Here we consider both pos- 
itive and negative impurity SPs. A single impurity is put at 
site (16, 16). We define the on-site magnetic order parameter. 
Mi = (-i)''\Yjfi('^ipi ~ nifji). This definition is suitable for 
the typical (n, 0) SDW order, consistent with previous experi- 
ments-'^ and theoretical calculations.Q2l2li 

The intensity plots of the site-dependent particle number 
= Zi//(«iA(T + "'Aii) ^'^'^ magnetic order M, in real space are 
shown in Fig. [3] The left panels of Fig. |3] plot the spatial 
distribution of the particle number. For positive SP, electrons 
are repelled by the impurity, therefore at the impurity site the 
value of particle number is reduced. Increasing the positive 
SP will lead to smaller values of n, at the impurity site, which, 
when Vs > 6, will vanish. For negative SP, on the contrary, 
electrons are attracted to the impurity and larger IV,! will lead 
to a higher particle number at the impurity site. The particle 
numbers will recover to the bulk value 2.0 at about 2 lattice 
constants away from the impurity site for both positive and 
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FIG. 3: (Color online) The intensity plots of the particle number (left 
panels) and magnetic order (right panels) at zero doping and zero 
temperature for different SPs Vs = 1,5,-1,-fOO. 



negative SPs. In doped samples, these characteristics of «, 
do not change except that the value of n,- far away from the 
impurity site will he 2 + x, where x is the electron doping 
concentration. 

The right panels of Fig. [3] show the real space modulation 
of the magnetic order M,. For small positive SP V, = 1, the 
values of magnetic order oscillate near the impurity site with 
the maximum M, = 0.1 15 at the impurity site, slightly higher 
than the bulk value 0.105. For moderate SP Vs - 5, the mag- 
nitude of Mi drops down almost to zero at the impurity site, 
with the maximum M, ~ 0.12 appearing in the vicinity of 
the impurity site. We thus expect that stronger SP will lead 
to stronger oscillation of M, around the impurity and this is 
proved by setting Vs = ±100 which we take Vs - -100 as an 
example as shown in Fig. [3] From the corresponding plot we 
can see that the modulation of M,- is very strong around the im- 
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FIG. 4: (Color online) The intensity plots of the LDOS at zero energy 
for different SPs V, = 3 and V, = 100. 
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purity, similar to Friedel oscillation and it will recover to the 
impurity-free value at about 4^5 lattice constants away from 
the impurity site. Order parameters for negative SP Vj = -1 
are also shown, unlike the enhanced M, in the - 1 case, 
we find that the magnitude of M, is suppressed at the impurity 
site and the oscillation of M, is stronger than the V, - 1 case. 

We now study the low-energy impurity-induced bound 
states. The intensity plots of the LDOS in real space at zero 
energy are shown in Fig. |4]for = 3 and = 100, re- 
spectively. As seen, for = 3 the LDOS at the impurity 
site (16, 16) is finite while for nearly unitary SP Vs - 100 
it vanishes. The LDOS modulates near the impurity site and 
some bright spots can be seen clearly around the impurity, in- 
dicating the existence of bound states at low energy. Another 
prominent feature revealed from the LDOS map in Fig. |4]is 
the four-fold symmetry breaking which is more obvious near 
the impurity. The symmetry of the system reduces to C2 and 
it survives for various SPs no matter they are positive or neg- 
ative. Furthermore, there also exists one-dimensional modu- 
lation of the LDOS along the y-axis even when it is far away 
from the impurity. This feature is similar to the experimen- 
tally observed nematic electronic structure,^^** thus supporting 
the impurity effect as a possible candidate for the formation 
of nematic order. Since in the two-orbital model we use, each 
unit cell contains two inequivalent Fe atoms, the existence of 
an impurity on either site of the unit cell will naturally break 
the four-fold symmetry of the system, thus we conclude that 
the breaking of the four-fold symmetry in the LDOS is in- 
duced not only by the SDW order, but also by the intrinsic 
asymmetry pinned by the impurity. 

Then we study the energy and site dependance of the 
LDOS. Without impurity, the LDOS is uniform and site inde- 
pendent, with its minimum located at negative energy, consis- 
tent with STM experiment.^'' There are two coherence peaks 
at negative energies and two at positive energies, as shown in 
Fig.|5] Since the LDOS has finite value at the impurity site for 
weak SP, we plot it exactly at the impurity site for both posi- 
tive and negative SPs in Fig.|6] For positive impurity potential 
Vs - 3 the spectrum displays two distinct resonance peaks at 
negative energies which are denoted by two arrows. The in- 
tensity of the left peak is higher than that of the right one. The 
splitting of the resonance peaks is due to the presence of the 
inter-orbital coupling 14 and the resonance peaks are related 
to the opening of SDW gap. As Vs is increased the two reso- 



FIG. 5: The bulk LDOS in the SDW state of undoped sample. 
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FIG. 6: (Color online) (a) The LDOS at the impurity site as a function 
of to, for different positive SPs V, = 3, 4, 5. (b) is similar to (a), but 
for negative SPs Vs = -3, -4, -6. (c) and (d) [(e) and (f)] are similar 
to (a) and (b), respectively, but are plotted at site (15, 17) [(15, 16)]. 



nance peaks shift to higher energies and the intensities of the 
peaks decrease as shown in Fig.|6](a). At last, the intensity of 
the right peak becomes higher than that of the left one. For 
Vs - 6, the LDOS at the impurity site nearly vanishes. The 
feature of the LDOS for negative SPs shown in Fig. |6](b) is 
different from that in the positive SP case, for example, the 
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intensities of the peaks are much lower. For Vj — -3, double 
peaks show up at both sides of the Fermi energy, which we 
also denote by two arrows. These peaks shift to lower ener- 
gies with increased value of IV,]. As V, reaches to Vs = -8, 
the LDOS at the impurity site also vanishes. 

In Figs. |6] (c) and |6] (d) we plot the LDOS at nnn site 
(15, 17) of the impurity for positive and negative SPs, respec- 
tively. One can see that with increased strength of SP, the dou- 
ble peaks move to higher (lower) energies for positive (neg- 
ative) SP. For positive SP, increasing impurity strength will 
lead to increased peak intensities and this is in contrast to that 
at the impurity site. However, for negative SP, the situation is 
similar to that at the impurity site. Figs.|6](e) and |6](f) plot 
the LDOS at nn site (15, 16). It is shown that the intensities 
of the impurity-induced resonance peaks are much lower than 
those at the impurity site (16, 16), although the characteristics 
are similar. Since the impurity has four nnn and nn sites, and 
the system has only C2 symmetry, there are two inequivalent 
nnn and nn sites, respectively. The LDOS at the other nnn 
site (15, 15) and nn site (16, 15) does not show the impurity- 
induced resonance peaks at low energies (not shown here) 
and resembles the bulk LDOS, again suggesting the four-fold 
symmetry breaking. 

The properties of the low-energy bound states shown in 
Fig. are significantly different from those in the pure SC 
state ™ In the pure SC state, the positions of the bound 
states are close to the SC coherence peaks so that they may be 
masked by the SC coherence peaks and may be hard to detect 
by experiments. However, in the SDW state, the energies of 
the bound states are close to the Fermi energy so that they can 
be easily detected by experiments. 



POSITIVE IMPURITY SCATTERING IN DOPED 
SAMPLES 
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FIG. 7: (Color online) The intensity plots of the SC [(a) and (c)] and 
magnetic [(b) and (d)] orders for weak SP V , = 1, at doping x = 0.04 
and X = 0.08, respectively. 
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As we have discussed above, in both pure SDW state and 
pure SC state, the bound states are induced by a single non- 
magnetic impurity. Since the detailed features of the reso- 
nance peaks are quite different between these two cases, the 
impurity effect in the underdoped regime where the SC and 
SDW orders coexist is an intriguing question. In particular, 
both theoretical analyses'^ and experimental observationfl^^ 
do suggest the coexistence of these two orders in this regime. 
In the following, we will not plot the real space particle num- 
ber since it is similar to the undoped case. We want to mention 
that around the moderate doping x = 0.08, the impurity could 
induce a weak charge density wave for various SPs. However, 
the 6n/{2 + x) is less than 0.5/100, so we neglect it. 

For small SP Vs - 1, we can see from Fig. [7] (a) that at 
low doping X - 0.04, the amplitude of the SC order A, is re- 
duced at and around the impurity site, which will recover to 
the impurity-free value at about 6 lattice constants away from 
the impurity. But the SC order is not always suppressed at 
the impurity site. As doping is increased to x > 0.08, at the 
impurity site the magnitude of the SC order is enhanced [see 
Fig.|7](c)], which means the impurity is not a pair breaker in 
this case. At jc = 0.04, the magnetic order M, at the impurity 



FIG. 8: (Color online) The LDOS at the impurity site as a function 
of aifor V, = 1,3,5, at different doping levels x = 0.04,0.08,0.12. 



site is enhanced, similar to that in the undoped case and we 
notice that there exist modulations along the diagonal direc- 
tions as can be seen from Fig.[7](b). At x = 0.08, the pattern 
of magnetic order changes, the system separates into two sub- 
lattices explicitly. The value of M,- in one sublattice is about 
0.05, while in the other one is ~ 0.007, with M, ~ 0.01 at the 
impurity site. This impuiity-induced two-sublattice pattern of 
magnetic order survives until the doping level is beyond the 
region where the SDW and SC orders coexist. At higher dop- 
ing X - 0.12, A, is enhanced just like the x = 0.08 case, but 
with a vanishingly small value of M,. 

For larger SP Vs = 3.0, the order parameters are similar to 
those for = 1 .0, except that M, is reduced at the impurity 
site at all doping levels. As doping increases to x - 0.08, the 
system also separates into two sublattices. An enhanced SC 
order will appear at the impurity site when x reaches to 0.12 
(not shown here). 

The LDOS has finite values at the impurity site for small SP 
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FIG. 9: (Color online) For V, = 1, the LDOS on nn and nnn sites of 
the impurity as a function of o) at various doping levels. The black 
dashed line represents the bulk LDOS. 



in all doped samples, but unlike the undoped one, the double 
resonance peaks are absent. For Vs - I, the effect of the im- 
purity on the LDOS is shown in the left panel of Fig. |8] In this 
case, at jc = 0.04, the intensities of the SC coherence peaks 
at both positive and negative energies are suppressed by the 
impurity. On the other hand, at doping x > 0.08, the negative 
SC coherence peak is enhanced by the impurity while the pos- 
itive one remains almost unchanged. For larger SP V, = 3, at 
X - 0.04, the intensities of the SC coherence peaks are further 
suppressed. When x = 0.08 there is a sharp in-gap resonance 
peak located at negative energy and close to the SC coher- 
ence peak [see Fig.[8]^b)]. As doping is increased to x = 0.12, 
the intensity of the in-gap peak becomes higher. On the other 
hand, for moderate SP Vs = 5, at x = 0.04, there exists an in- 
gap bound state at negative energy while at both x = 0.08 and 
X - 0.12, there are two in-gap bound states, one at positive 
energy, the other one at negative energy. The magnitude of 
the LDOS at all doping levels becomes considerably smaller 
and reaches to zero for larger SP. 

For Vs = 1 .0, the impurity induces only minor modulations 
on the LDOS around the impurity site, which is similar to 
the bulk LDOS at all doping levels [see Fig. |9]. The positive 
energy peak at nn and nnn sites is higher than the negative 
one at low doping x - 0.04. As x increases, the intensity of 
the negative resonance peak gradually becomes higher than 
that of the positive one, similar to the evolution of the bulk 
LDOS with doping. Although the system does not have C4 
symmetry, the main features of the LDOS at the four nnn (nn) 
sites are similar to each other, thus in Fig.|9j we only plot the 
LDOS at one of the nnn (nn) sites for clarity. 

As the SP increases to Vs - 3.0, we show the LDOS at 
two inequivalent nnn (nn) sites in the upper (lower) panel of 
At low doping x - 0.04, the effect of the impurity is 
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Fig. 10 



weak and no in-gap bound states exist at nnn and nn sites. At 
higher doping x = 0.08, a single in-gap resonance peak shows 
up at both the two nnn sites with different intensities, but their 
positions are similar to each other, both are located at positive 
energy and close to one of the SC coherence peaks. There 
is also a single in-gap peak at the two nn sites, the LDOS 
of which is identical to each other As doping increases to 
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FIG. 10: (Color online) For Ki = 3, the LDOS on nn and nnn sites of 
the impurity as a function of ai at various doping levels. The upper 
panel is for the two nnn sites (15, 15) (red solid) and (15, 17) (pink 
short dot) and the lower one is for the two nn sites (15, 16) (green 
solid) and (16, 15) (blue short dot). The black dashed line represents 
the bulk LDOS. 



X - 0.12, the LDOS at the nnn sites is similar to the x - 0.08 
case, except for a higher peak intensity at positive energy and 
the addition of a hump at negative energy on one of the nnn 
sites, which will evolve into a resonance peak when further 
increasing doping (not shown here). The LDOS at the two 
nn sites is also identical to each other and clearly shows two 
in-gap resonance peaks. 

For moderate SP = 5, the difference of the LDOS be- 
tween the two inequivalent nnn sites (15, 15) and (15, 17) be- 
comes remarkable [see the upper panel of Fig. [TT|. At one 
of the nnn sites the single in-gap resonance peak is located at 
positive energy while at the other one it is located at negative 
energy. As doping increases, the peaks at both the two nnn 
sites move closer to the Fermi energy. On the other hand, at 
the two nn sites (15, 16) and (16, 15), there is a single in-gap 
resonance peak located above the Fermi energy at all dopings. 
At low doping the LDOS spectra at the two nn site are differ- 
ent. As doping increases the peaks shift to the Fermi energy 
and the LDOS at the two nn sites will become identical. 

For nearly unitary positive SP = 100, the SC and mag- 
netic orders are both suppressed and oscillate in the vicinity 
of the impurity site and their magnitudes reach the minimum 
exactly at the impurity site. The suppressed order parameters 
recover to their bulk value at about 3 lattice constants away 
from the impurity. We note that when x > 0.04 the system 



8 



1.6 



0.8 



0.0 



x=0.04 



x=0.08 



nnnl 


(a) 




bulkj 








Hi (1 



















(b) 












f 



3 



1.6 



0.8 



0.0 



■0.2 0.0 0.2 

(0 



x=0.04 



-0.2 0.0 0.2 

CO 



x=0.08 




(d) 



f : I 
/I in. : 
/J i;i ■■ 



.V 




■0.2 0.0 0.2 

CO 



-0.2 0.0 0.2 -0.2 0.0 0.2 



CO 



CO 



FIG. 11: (Color online) Similar to Fig.fTO] bur for V, = 5. 
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FIG. 12: (Color online) Similar to Fig.|7] but for V, = 100. 



will separate into two sublattices. The difference of M, be- 
tween the two sublattices is larger at x - 0.08 than that at 
X = 0.04 while the magnitude of M, decreases with doping 
and will vanish at x > 0. 1 . 

Since the LDOS at the impurity site is zero for such a strong 
SP, we thus plot it on nnn and nn sites of the impurity in 
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FIG. 13: (Color online) The LDOS on nn and nnn sites of the im- 
purity as a function of a>, for Vs = 100. Different panels correspond 
to different doping levels. The red solid line denotes the nnn site 
(15, 17) while the blue dashed line denotes the nn site (16, 15), with 
the bulk LDOS denoted by the black dashed line. 



resonance peak appears close to zero energy on the positive 
side. It shifts to negative energy with reduced intensity as x 
increases to 0.04. When x > 0.05 two in-gap resonance peaks 
show up. As doping increases further, they are pushed away 
by each other and finally merge into the SC coherence peaks 
of the bulk LDOS at x > 0.08. In all cases, the LDOS exhibits 
clear C2 symmetry. However, the LDOS on the inequivalent 
nn (nnn) sites is qualitatively the same, thus we choose the 
sites (16, 15) and (15, 17) as an example for convenience. 



VI. NEGATIVE IMPURITY SCATTERING IN DOPED 
SAMPLES 

In real materials, both positive and negative SPs are possi- 
ble, and the response of the system to the impurity may de- 
pend on the sign of the SP, thus we discuss the negative SP 
case in this section. 

As we can see from Fig. 



14 for y, 



-3, the impurity 
induces oscillation of the SC order with A, being suppressed 
at the impurity site and enhanced on several nearby sites at x = 
0.04. At higher doping x = 0.08, around the impurity site A, is 
suppressed and the oscillation is not distinct. The magnitude 
of Mi is suppressed at the impurity site at all doping levels 
and apparently M, will divide into two sublattices at relatively 
higher doping. At x - 0.08, at the impurity site M, is close to 
zero, and is 0.05 and 0.01 on the two sublattices, respectively. 

Contrary to the = 3 case where the impurity effect is 
stronger at relatively higher doping, the in-gap resonance peak 
for negative SP is clearer at lower doping. From Fig. 15 we 



Fig 13 at various dopings. At x = 0.02, a sharp in-gap can see, for Vs = -3, at the impurity site there is a sharp 
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FIG. 14: (Color online) Similar to Fig.|7] but for V, = -3. 
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FIG. 15: (Color online) The LDOS at the impurity site as a function 
of oj at different dopings x = 0.04, 0.08, 0. 12, for V, = -3, -5, -8. 



in-gap resonance peak at positive energy at low doping x - 
0.04. At higher doping, the LDOS at the impurity site has two 
peaks at the edge of the SC coherence peaks with the right 
peak being higher than the left one. Here we do not show 

we 



the SC coherence peaks in Fig. 15 but in Figs. 16 anc 18 
plot the corresponding bulk LDOS in which the SC coherence 
peaks are explicitly shown. As the SP strength increases to 
Vs - -5, the corresponding peaks are all suppressed and in 
the X - 0.04 case, the in-gap resonance peak shifts towards 
the Fermi energy. From the Fig. 15 b) we can see that the 
intensity of the LDOS at the impurity site nearly vanishes for 

y, = -8 . 

Then we plot the LDOS on nnn and nn sites for Vs - -3 in 
Fig. [16] at different dopings. At x = 0.04, the LDOS on both 
the nnn and nn sites shows the existence of in-gap resonance 
peaks, which gradually merge into the SC coherence peaks as 
doping increases. Again, the breaking of four-fold symmetry 
is more obvious on the nnn sites than it is on the nn sites and 
as doping increases, this asymmetry tends to diminish. 
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FIG. 16: (Color online) For K, = -3, the LDOS on nnn and nn sites 
of the impurity as a function of a) at various dopings. The upper and 
lower panels are for nnn and nn sites, respectively. The bulk LDOS 
is denoted by the black dashed line. 



For strong negative potential Vj — -8, both the SC and 
magnetic orders are suppressed at the impurity site and oscil- 
late around it [see Fig.[T7). At about 4 lattice constants away 
from the impurity, the SC and magnetic orders recover to their 
bulk values. At doping x - 0.08, M, also separates into two 
sublattices as can be seen from Fig.[T7](d). 



As shown in Fig. 18 



there are in-gap resonance peaks on 
all the nnn and nn sites at x = 0.04. The breaking of the four- 
fold symmetry is minor, although still visible. The intensities 
of the in-gap peaks are suppressed with increased doping. At 
X - 0.12 no in-gap peaks exist. The profile of the LDOS on 
nn sites is similar to that on nnn sites, except that the coiTe- 
sponding peaks are lower and the asymmetry is weaker. 

For nearly unitary negative potential - -100, the real 
space distributions of the SC and magnetic orders are similar 



to the positive unitary potential case. Fig. 19 shows the LDOS 
on nnn and nn sites at various dopings. Similar to the strong 
positive SP case, there is a sharp in-gap peak at low doping 
X - 0.04, the position of which is almost at the Fermi energy. 
We believe that this sharp peak is due to the existence of the 
SDW order. As doping increases to x = 0.08, the height of 
the in-gap peak drops and becomes lower than that of the SC 
coherence peak. As doping increases to x = 0.12, the in-gap 
peak disappears with two resonance peaks being very close to 
the SC coherence peaks. 
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FIG. 17: (Color online) Similar to Fig.|7] but for V, = -8. 
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FIG. 18: (Color online) Similar to Fig.fTel but for V, = -8. 



VII. SUMMARY 

By solving the BdG equations self-consistently, it is shown 
that without impurity, at zero doping and in the SDW state, 
there exist equal-sized electron- and hole-like FS pockets 
along the F - M line of the BZ, inside which the Dirac cones 
form. The electron- and hole-like Dirac cones appear in-pairs 



near the Fermi energy and are located very close to each other. 
The effect of doping is mainly to reduce the size of the hole 
pockets while increase that of the electron ones, consistent 
with the increased electron density. 

When impurity is introduced into the system, we find that in 
the parent compound, strong SP, no matter repulsive or attrac- 
tive, could induce considerably large oscillation of the mag- 
netic order around the impurity site. In addition, for all the SP 
strength we investigated, there exists one-dimensional mod- 
ulation of the LDOS, similar to the experimentally observed 
nematic electronic structure, thus supporting the impurity ef- 
fect as a possible candidate for the formation of nematic order 
Furthermore, two impurity-induced resonance peaks are ob- 
served around the impurity site and they are shifted to higher 
(lower) energies as the strength of the positive (negative) SP 
is increased. 

In doped samples, generally speaking, the SC and magnetic 
orders are suppressed at and around the impurity site, with 
more complicated variations compared to those in cuprates. 
However, for positive SP at higher doping or negative SP 
at lower doping, the SC order may even be enhanced at or 
around the impurity site, suggesting that the impurity is not 
a pair breaker in this case. In addition, impurity could sepa- 
rate the system into two sublattices denoted by two different 
values of magnetic order which can be seen more clearly at 
relatively higher doping. Furthermore, there exist impurity- 
induced bound states at and around the impurity site, whose 
positions and numbers depend on the strength and sign of the 
SP, as well as on the doping concentration. For weak and 
moderate SPs, a distinct bound state exists explicitly at the 
nnn sites of the impurity. For unitary impurity, there is a sharp 
in-gap peak at low doping, while at high doping, the impurity 
induced bound state is close to the SC coherence peaks. On 
the other hand, in a small range of moderate doping there are 
two in-gap peaks only for positive SP. 

In all cases, impurity breaks the four-fold symmetry of the 
system and has a stronger effect on nnn sites than it does on nn 
sites as can be seen from the LDOS. This symmetry breaking 
is induced not only by the SDW order, but also by the intrinsic 
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asymmetry in our model pinned by the impurity. All the above 
features could be used to detect the presence of the SDW order 
and to probe the coexistence of the SDW and SC orders. 
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